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SUMMARY 

In this study we address the question under which conditions a saturated velocity field stem- 
ming from geodynamo simulations leads to an exponential growth of the magnetic field in a 
corresponding kinematic calculation. We perform global self-consistent geodynamo simula- 
tions and calculate the evolution of a kinematically advanced tracer field. The self-consistent 
velocity field enters the induction equation in each time step, but th e tracer field does not con- 
tribute to the Lorentz force. This experiment has been performed by Catt aneo & Tobias] d2009l) 
and is closely related to the test field method by Schrinner et al. (2005, 2007). We find two dy- 
namo regimes in which the tracer field either grows exponentially or approaches a state aligned 
with the actual self-consistent magnetic field after an initial transition period. Both regimes can 
be distinguished by t he Rossby number and coinc ide with the dipolar and multipolar dynamo 
regimes identified bv lChristensen & Aub ert (2006). Dipolar dynamos with low Rossby number 
are kinematically stable whereas the tracer field grows exponentially in the multipolar dynamo 
regime. This difference in the saturation process for dynamos in both regimes comes along 
with differences in their time variability. Within our sample of 20 models, solely kinematically 
unstable dynamos show dipole reversals and large excursions. The complicated time behaviour 
of these dynamos presumably relates to the alternating growth of several competing dynamo 
modes. On the other hand, dynamos in the low Rossby number regime exhibit a rather simple 
time dependence and their saturation merely results in a fluctuation of the fundamental dynamo 
mode about its critical state. 

Key words: Dynamo: theories and simulations; Earth's core; geomagnetic field; magnetohy- 
drodynamics. 



lation. Despite the fact that the magnetic field is saturated in the 
full non-linear system, it can grow in a kinematic treatment, be- 
cause both associated linearized stability problems are different. 
The flows taken from a saturated dynamo simulation and then used 
in a kinematic calculation need only quench the growth of the par- 
ticular magnetic field found in the nonlinear problem and can in 
principle allow others to grow. As lTilgner~ & Brandenburg (2008) 
have pointed out there is at least one example, the benchmark dy- 
namo case 1 dChristensen et alj2001f) . where the field taken from a 
saturated dynamo is also kinematically stable. 

In this study, we show that there is in fact a whole class of satu- 
rated, chaotic, time-dependent dynamos for which the correspond- 
ing kinematic dynamo is stable. In order to assess kinematic stabil- 
ity - in the sense explained above - we solve the MHD-equations 
for a Boussinesq fluid in a rotating spherical shell. At the same time 
we evolve a second passive tracer field using the induction equa- 
tion. While the tracer field experiences the self-consistent velocity 
field at each time step it does not contribute to the Lorentz force. 
This method has been used bv lCattaneo & Tobias! J2009I) for box 
simulations and a shell model and is closely related to the test-field 



1 INTRODUCTION 

The time variability of cosmic magnetic fields has always been an 
argument in favour of hydromagnetic dynamo action. Its under- 
standing is crucial for insights in the interior dynamics of stars and 
planets. The time dependence of convective dynamos is attributable 
to a non-stationary buoyancy flux as well as to a time dependent 
equilibration of the magnetic field. The latter is subject of the study 
presented here. 

How do dynamos saturate and in particular in which way is the 
saturation reflected in their time dependence? In a general descrip- 
tion, the infinite growth of a magnetic field due to an appropriate 
motion of a conducting fluid is inhibited owing to the backreaction 
of the Lorentz force on the flow; the resulting changes in the flow 
cause a reduction of dynamo action. Flows which are influenced 
by t he Lorentz force in this w ay are call ed saturated. Neverthe- 
less. ICattaneo & Tobiasl (2009) as well as |jilgner & Brandenburg! 



(2008) demonstrate that saturated flows may lead to exponential 
growth of the magnetic field in a corresponding kinematic calcu- 
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metho d to determine mean-field coefficients dSchrinner et alj2005l 
l2007b . 

Within a sample of 20 models, we identify two dis- 
tinct dynamo_j«gimes_d£pendent on a modified Rossby num- 
ber dChristensen & Aubertll200^ in which the tracer field either 
grows exponentially or reaches a state aligned with the actual self- 
consistent magnetic field after an initial transition period. More- 
over, differences in the kinematic stability of the dynamos are 
linked to differences in their time variability: Exclusively kinemat- 
ically unstable dynamos in the high Rossby number regime show 
polarity reversals of the axial dipole field. We attribute the com- 
plicated time behaviour of these models to an alternating growth 
of many competing dynamo modes. On the other hand, the eigen- 
value computation suggests that dynamos with low Rossby number 
are dominated by only one fundamental mode which is repeatedly 
quenched and rebuilt. All other modes in this case are clearly sub- 
critical. In this sense, dynamo models in the low Rossby number 
regime, i.e. fast rotators, exhibit a simple time dependence and their 
time- variability consists of fluctuations about their critical state. 



2 DYNAMO CALCULATIONS 

We consider an electrically conducting Boussinesq fluid in a ro- 
tating spherical s hell and solve the MHD-e quations as given by 
l01sonetal.l ( fl999h and described in detail bv lChristensen & Wichtl 
(2007). In addition, we compute the evolution of a passive tracer 
field with the help of a second induction equation 



dB Tl /dt = V x (u x Bti) + l/Pm\7 2 B 1 



(1) 



While the tracer field, Bti, experiences the self-consistent velocity 
field u in each time step, it does not contribute to the Lorentz force. 
Hence it does not act on the velocity field and is "passive" in this 
sense. The initial conditions for the tracer field have been chosen 
arbitrarily with the help of a random number generator. Moreover, 
for models 10-15, we added some random noise to the tracer field 
in each time step. This enables us to perturb the tracer field perma- 
nently and prevents it from becoming aligned with the actual, self- 
consistent field. In these simulations, we advance the tracer field for 
at least 10 magnetic diffusion times in order to test for kinematic 
stability. 

According to the scaling we used, the equations are governed 
by four parameters. These are the Ekman number E = v /Q,D 2 , the 
(modified) Rayleigh number Ra — argoATD / Vf2, the Prandtl 
number Pr = vjn and the magnetic Prandtl number Pm = v/rj. 
In these expressions, u denotes the kinematic viscosity, Q. the ro- 
tation rate, D the shell width, «t the thermal expansion coeffi- 
cient, go is the gravitational acceleration at the outer boundary, AT 
stands for the temperature difference between the inner and outer 
spherical boundaries, k is the thermal and r\ — 1/ficr the mag- 
netic diffusivity with the magnetic permeability ji and the electrical 
conductivity a. All four parameters have been varied to build up a 
sample of 20 dynamo models, see Table 1 . 

Output parameters used here in order to interpret the results 
are the magnetic Reynolds number, Rm = UD/r/, the Elsasser 
number, A = B 2 /q(at]£1, and the Rossby number, Ro — U/Dfl. 
In these expressions, U and B denote rms-values of the velocity 
and the magnetic field inside the shell, respectively, and p is the 
density. Furthermor e, we adopt the defini t ion of a local Rossby 
number proposed bv lChristensen & Auberj J2006l) . 
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Figure 1. Magnetic energy densities for two computational runs of model 2. 
Both runs have been started from very similar initial conditions which differ 
only by a small deflection (dashed line) of the magnetic dipole axis. Nev- 
ertheless, both models evolve differently which demonstrates the chaotic 
character of these dynamos. 



Here, l/nis the mean half wavelength of the flow and I is the mean 
harmonic degree derived from the kinetic energy spectrum, 

, < ui ■ m > 



< u ■ u > 



(3) 



Roi = Ro ■ — 



(2) 



The brackets in Eq. (O denote an average over time and radii, u% 
stands for the velocity component of harmonic degree I. 



3 RESULTS 

Within our 20 examples (see Table [TJ we find 5 dynamos which 
are kinematically unstable and 14 which are kinematically stable. 
One example (model 15) belongs to both classes; although in gen- 
eral unstable, the tracer field does not grow within certain periods 
of several magnetic diffusion times. Note that all dynamos con- 
sidered here operate in the so called strong field regime, i.e. the 
Elsasser number is of order unity or larger. The equatorial symme- 
try is broken for most of the kinematically stable and all unstable 
models. Except for mod el 1, the quasi-steady benchmark dynamo 
dChristensen et alj200ll) . all models exhibit highly time-dependent 
or even chaotic behaviour. This is demonstrated in Fig.[T]for model 
2, the next simplest example to the benchmark dynamo. This dy- 
namo appears to be chaotic, and as an experiment we performed 
two simulations starting from almost identical initial conditions 
(the difference between two initial conditions is a small deflection 
of the magnetic dipole axis in the second run). The evolution from 
both initial conditions is shown in Fig.[T]where the magnetic energy 
densities can be seen to diverge rapidly with time. 

The regimes of kinematically stable and unstable dynamos can 
be clearly distinguished by the modified Rossby number (see Table 
[]}, Roi. Models with low Rossby number are kinematically sta- 
ble whereas the tracer field grows exponentially for dynamos in the 
high Rossby number regime. The transition between both regimes 
occurs at Roi ~ 0.12. There are two further properties related to 
both regimes which deserve mentioning. All dynamos we found 
to be kinematically stable are dipolar and do not show any polar- 
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Figure 2. Relative dipole field strength /di p versus local Rossby number 
Roi . Stars denote non-reversing dynamos which are kinematically stable, 
whereas triangles represent dynamos which do reverse and are kinemati- 
cally unstable. Both regimes coincide with the dipolar and multipolar dy- 
namo regimes identified by Christensen and Aubert (2006). There is one 
example (diamonds), model 15, which undergoes a transition between both 
regimes. Note that this example has a considerably lower relative dipole 
field strength in its second state. 



ity reversals, while dynamos in the second regime are multipolar 
and do reverse. This is also illustrated in Fig. [2] Here, the relative 
dipole field strength, /dip, on the outer shell boundary is plotted 
versus the modified Rossby number, Roi\ /dip is defined as the 
time-average ratio of the dipole field strength to the field strength 
inharmonic degrees 1 to 12. Bo th regimes visible in Fi g.[2]coincide 
with those identified earlier bv lChristensen & Aubertl d2006h . Fig- 
ure[3]compares the magnetic energy densities of the tracer field for 
a kinematically stable (model 8) and a kinematically unstable dy- 
namo (model 19), varying with time. While the tracer field grows 
rapidly after an initial transient phase in the latter case, it reaches 
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Figure 4. Snapshot of the radial component of the actual magnetic field 
and the tracer field for model 8 taken some time after an initial transition 
period at r = 0.62r o where r is the outer shell radius. Note that the tracer 
field is completely aligned with the actual magnetic field. Both components 
are normalised due to their maxima and minima. Therefore the greyscale 
coding varies from -1, white, to +1, black, and the contour lines correspond 
to ±0.1, ±0.3, ±0.5, ±0.7, ±0.9. Following contour plots are presented 
in the same style. 



a state aligned with the actual field if the dynamo is kinematically 
stable. Then, the energy density of the tracer field normalised with 
the energy density of the actual self-consistent field approaches a 
constant level which depends only on the initial conditions. This 
is also confirmed by looking at the corresponding field configura- 
tions. Figure[4]displays the radial component of the tracer field for 
model 8, which differs from the actual field only by an overall scale 
factor. Therefore, only one contour plot is given. On the other hand, 
although they have similar spatial scales, both field components are 
clearly not aligned but very different for model 19 (see Fig.[5}- 



Radial Magnetic Field 




Radial Tracer Field 




Figure 3. Ratio of the magnetic energy densities for the tracer field, 
i? magTr , and the actual magnetic field, -E mag , versus time for model 8 
(dashed line) and model 19 (solid line). 



Figure 5. Snapshots of the radial component of the actual magnetic field 
(top) and the tracer field (bottom) for model 19 at r = 0.62r o . Contour 
lines: see figure|4] 
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Table 1. Overview of the runs considered, ordered with respect to their modified Rossby number. All kinematically unstable 
models exhibit dipole reversals whereas all kinematically stable models do not. 
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unstable state (see Fig. [6]l and vice versa. As long as the tracer 
field remains stable, the tilt angle of the dipole axis fluctuates about 
the actual polarity state. However, when the tracer field becomes 
unstable, also the polarity of the dipole field starts to reverse. This 
coincidence is observed for transitions in both directions, i.e. the 
tilt angle of the dipole axis also stabilises when model 15 becomes 
intermittently stable. While the magnetic field is quite dipolar with 
/dip = 0.61 for periods in which the polarity and the tracer field 
are stable the relative dipole field strength decreases drastically to 
/dip = 0.16 otherwise. The strong connection of field morphology, 
time dependence and saturation is not only present separately in 
several models but manifests itself in the time variation of a single 
dynamo model, too. 



4 DISCUSSION 

The existence of kinematically unstable dynamos wa s expected 
dCattaneo & Tobiasll2009l ; lTilgner & Brandenburdl2008l) . The find- 
ing of a class of kinematically stable but yet time-dependent or 
even chaotic dynamos, however, needs some further explanation. 
The lack of growing modes for these models already suggests that 
almost all field configurations for the tracer field are decaying, ex- 
cept the one aligned with the actual, self-consistent field. But this 
component of the tracer field is quenched by the saturated veloc- 
ity field. Thus, the tracer field follows the actual field with time, 
apart from a different, arbitrary amplitude due to the linearity of 
the induction equation. 

This interpretation is confirmed by looking at the spectrum of 
the time and azimuthally averaged dynamo operator D, 

Db l = \ l b l (4) 
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Figure 6. Tilt angle of the dipole axis for model 15 as a function of time 
(solid line) and magnetic energy density of the tracer field normalised by 
the magnetic energy density of the actual field, E maETr /E mag (dashed 
line which runs out of the figure at roughly 2.3 magnetic diffusion times). 
As soon as the dynamo reverses it becomes kinematically unstable. 



Model 15 is in general kinematically unstable but also exhibits 
periods of several magnetic diffusion times in which the tracer field 
stays stable. According to its local Rossby number, Roi — 0.114, it 
is located close to the boundry between both dynamo regimes and 
undergoes transitions from one to the other. 

We could detect transitions from a kinematically stable to an 
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Figure 7. Radial components of the thirst three dipolar eigenmodes b l ,i = 
1 ... 3 of the time averaged dynamo operator for model 2. The correspond- 
ing eigenvalues are Ai = —3.87, A2 = —34.83 and A3 = —42.45 in 
units of r)/D 2 . Note the huge drop in decay rates after \±. 
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with eigenmodes b % and eigenvalues A ! . In this, the operator D is 
defined as 

Db = V x (u x b + a ■ b - (3 Vb - 77V x b) (5) 

Note that D, also known as mean-field dynamo operator 
( lKrause&RadleJl980l) . contains the mean velocity field u as well 
as the so called mean-field coefficients a. and /3, w hich are tensor s 
of second and third rank, respectively. As noted by Hovng (2009), 
these quantities appear inevitably as a consequence of averaging. 
They depend on the velocity field and the magnetic diffusivity of 
the considered dynamo model onl y and have been determined w ith 
the help of the test field method l lSchrinner et al. I I2005L |2007|) . A 
detailed discussion on the applicability of mean-field concepts to 
direct numerical simulations of rotating magnetoconvection an d a 
(quasi-)stationary dynamo is provided by Sc hrinner et alj j2007l> . A 
similar discussion for time dependent dynamos is not given here but 
will be subject of a forthcoming paper. A recent revie w on the test- 
field m ethod and its applications has been given by iBrandenburd 
d2009t> . 

Eigenv alues and eigenmodes of D have been computed as 
reported by 1 Schrinner et alj d2009h for model 2. In Fig. [7] the ra- 
dial components of the first three (dipolar) eigenmodes, b 1 , are dis- 
played. All modes decay exponentially; this had to be expe cted for 
kinem atically stable dynamos (see also the discussion in Hovng 
2009). However, the decay rates are given here in units of r//D , in 
which the molecular diffusivity r\ is about 30 times smaller than the 
turbulent one inferred from components of f3. Thus, l/|Ai | w 1/4 
is much larger than one effective diffusion time and the first, funda- 
mental, eigenmode is indeed close to its critical state. Due to a not- 
icable gap in decay rates after the fundamental mode, | Ai j <C | A2 1 , 
this is not equally true for the subsequent eigenmodes. They are 
much more diffusive, thus leaving the fundamental mode as the 
preferred eigenstate of the dynamo. Hence, the time dependence of 
model 2 may be understood in parts as a fluctuation of the funda- 
mental mode about its critical state. 

The dominance of the first eigenstate is also revealed by a de- 
composition of the actual, time-dependent magnetic field of model 
2 in a set of eigenmodes b 1 of D, 

B(r,t) = JV^r) ^ 

i 

The time-dependent and in general complex mode coefficients 
a ! (t) have been computed as 



Figure 8. Axisymmetric magnetic energy density (solid line) and energy 
contribution of the first, fundamental eigenmode (dashed line) varying with 
time. The fundamental eigenmode a 1 (t)b 1 (r) already contributes 75% up 
to 85% to the total amount and its time variability reflects much of the time 
dependence of the axisymmetric magnetic field. 

«'(*)= I j\r)-A{r,t)&\ (7) 
Jv 

in which j denotes the adjoint of the current j l — V x b l , 
and A is the vector potential of the actual, time-dependent field, 
B = V x A. The integration is carried out ov er the whole fluid 
domain V . For a deriv ation of Eq. ((7) we refer to Hovng ( 200jJ and 
ISchrinneretaI] ( l2009l) . 

In Fig. [8] the energy contribution of the fundamental eigen- 
mode a 1 (i)b 1 (r) is compared with the total axisymmetric mag- 
netic energy density. The fundamental eigenmode contributes at 
least 75% up to 85% to the total amount, revealing again its per- 
manent dominance throughout the simulation. 

T he equilibration pro cess for model 2 has been studied ear- 
lier by lOlson et alj J 19991) . too. They found that in regions with 
high magnetic energy density, the Lorentz force simply reduces 
locally the flow velocity without changing the overall pattern of 
convection. They investigated possible changes in the velocity, if 
the magnetic field and thus the Lorentz force is arbitrarily reduced 
at some instant in time and then recovers towards its equilibra- 
tion value. The kinematic effects relevant for dynamo action identi- 
fied by them, an a-effect from helicity in the columnar convection 
and an anti-cj effect from the mean azimuthal flow, were present 
in the same proportions close and far from equilibrium conditions 
of the magnetic field. Their finding is supported by the study pre- 
sented here. Saturation may reduce the amplitudes of a and thus 
the growth rates of the eigenstates of the related dynamo operator, 
but does not change their relative order. Therefore, the preferred 
eigenstate stays the same throughout the simulation. This is clear 
as there is a large gap between the growth rate of the fundamental 
eigenmode of the time-averaged dynamo operator D and all other 
eigenmodes, as mentioned above. 

So far we only analysed model 2 in detail. Here the veloc- 
ity field is nearly symmetric with respect to the equatorial plane 
and the magnetic field belongs to the dipolar family. Contributions 
of quadrupolar type are not present. In a more complicated exam- 
ple with broken equatorial symmetry, we expect the fundamental 
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quadrupolar mode to be excited. Although its growth rate will be 
smaller than the one for the fundamental dipolar mode and prob- 
ably subcritical, it is typically of the same order. In such a case 
a clear dominance of only one fundamental dipolar mode can no 
longer be deduced from the spectrum of the time-averaged dynamo 
operator D, and a second, quadrupolar mode may become impor- 
tant. 

F or the models in th e high Rossby number regime the find- 
ings of iCattaneo & Tob ias ( 2009) apply. These models act kinemat- 
ically as dynamos and the dynamo operator D possesses in general 
growing eigenmodes. A kinematic treatment of these dynamos does 
not reveal their actual time dependence. However, the results pre- 
se nted here suggest that the reg ime of dipolar dynamos identified 
by IChristensen & Aubertl ( 120061) is kinematically stable. For these 
models, the quenching of any magnetic field is fully captured in 
the velocity field and a kinematic treatment may indeed reproduce 
their actual time-dependence. Models of this dynamical regime are 
applicable to plan etary dynamos and proba bly also to dynamos of 
fast rotating stars ICh ristensen et al. 2009a), thus covering a large 
range of magnetic Reynolds numbers. Hence, an attempt to explain 
the kinematic stability of these models due to a magnetic Reynolds 
number which is close to its critical value fails. We emphasise again 
that the transition between both regimes is governed by the local 
Rossby number Roi and not by Rm, as can be already seen from 
TableQ] In the low Rossby number regime, the rotational constraint 
leads to columnar structured flows, dipolar magnetic fields and fi- 
nally to a rather simple time dependence, although these models op- 
erate in general far away from the dynamo threshold at Rm c « 40. 

Dipolar dynamo models which show occasionally reversals 
are located close to the regime boundary in Fig. [2] with Roi < 
0.12. They resemble the geodynamo in many respects and are 
therefore of particular interest. Explai ning polarity reversals of an 
otherwise predominantly dipolar field, Olson & Christ ensenl d2006l) 
suggest that the geodynamo crosses the boundary towards the mul- 
tipolar dynamo regime from time to time. With the help of scal- 
ing laws derived from numerical models, they indeed succeed in 
predicting a local Rossby number of Roi ~ 0.09 for the Earth's 
core. Adopting this viewpoint we link the occurence of geomag- 
netic reversals to a change in the saturation process. The quenching 
of a previously dipolar field may result in the preference of differ- 
ent, higher order modes if inertia gains importance in comparison 
to the coriolis force, and the dynamo undergoes an excursion into 
the kinematically unstable regime. Subsequently the dipole field is 
built up again, but it may have either polarity. A computation of 
eigenmodes and a mode decomposition similar to ® for model 15 
seems to be a promising approach to confirm this picture. Note that 
from the viewpoint we take here, the existence of dipolar, stable 
periods for model 15 demands more explanation than the fact that 
it reverses. 

However, whether inertia is indeed as important for the geo- 
dynamo as it is for presen t dynamo models is under debate (e.g. 
ISreenivasan & Jo nes 2006). In fact, the assumption of Roi » 0.09 
for the Earth's core leads to a characteristic length scale of only a 
few hundr ed meters, on which the magnetic field would be highly 
diffusive ( IChristensen et alj|2009bl) . 



5 CONCLUSIONS 

Fast rotating dynamos, characterised by a low Rossby number, 
are kinematically stable. Within this regime, a saturated velocity 
field taken from dynamo simulations does not lead to exponential 



growth of the magnetic field in a corresponding kinematic calcula- 
tion. Hence, saturation may be understood as a quality of the ve- 
locity field, only. For these dynamos, saturation results in the un- 
changed preference of a fundamental eigenstate, whereas different 
eigenmodes may supersede each other if inertia gains importance. 
This difference in the saturation process involves differences in the 
morphology of the magnetic field and its time dependence. Kine- 
matically stable dynamos are dipolar and exhibit a rather simple 
time variability, which may be interpreted as the fluctuation of the 
fundamental mode about its critical state. Kinematically unstable 
dynamos are much more complicated. The alternating growth of 
various modes leads to a multipolar field morphology and polarity 
reversals of the dipole field appear as a natural consequence. 
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